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In cosmological models with a varying gravitational constant, it is not clear whether primordial 
black holes preserve the value of G at their formation epoch. We investigate this question by using 
the Tolman-Bondi model to study the evolution of a background scalar field when a black hole 
forms from the collapse of dust in a flat Friedmann universe. Providing the back reaction of the 
scalar held on the metric can be neglected, we hud that the value of the scalar held at the event 
horizon very quickly assumes the background cosmological value. This suggests that there is very 
little gravitational memory. 

PACS numbers: 04.70.Bw, 04.50.+h, 97.60.Lf, 98.80.Cq 


I. INTRODUCTION 

Scalar-tensor (ST) theories of gravity provide a natural alternative to general relativity (GR). They describe gravity 
with not only a metric gab but also a scalar field (j). Derivatives of cj) appear as source terms in the field equations 
and (j) itself satisfies a wave equation. The strength of the gravitational coupling is determined by the function 
where GR is recovered in the limits w —> oo and {duj/dcfi) 0. 

ST theories can also be regarded as being equivalent to GR with a varying gravitational “constant” G. The most 
simple example of such a ST theory is Brans-Dicke theory |^, where ui{(j)) is constant and G oc (j)~^. However, weak 
field experiments have shown that lu > 500 and so the deviation from GR is small. For more general ST theories, 
where w is not constant, it is possible that to was much smaller at earlier times. So observations allow such theories 
to greatly deviate from GR in the early universe. 

There has been a renewed interest in ST theories in recent years due to the effective low energy actions of string 
theory involving one or more scalar fields. These scalar fields enter the field equations in much the same way as the 
scalar field in ST theories |3|. Also, the increasing popularity of inflation and quintessence suggests that scalar fields 
might need to be incorporated into cosmological models. 

The purpose of this paper is to study the effect of an evolving scalar field on the formation and evolution of a 
primordial black hole. In an asymptotically flat spacetime it is well known that a black hole radiates away any 
inhomogeneities in the scalar field until it becomes a stationary solution with constant (j) Q. This is a consequence 
of the famous “no hair” theorem. However, in ST cosmological models the scalar field is evolving with time and this 
would modify how the black hole evolves during its lifetime. 

Barrow was the first to examine this problem. He considered the two extreme possibilities: scenario (A), where 
the scalar field evolves everywhere homogeneously in the same way as the cosmological background; and scenario (B), 
where the black hole forces the scalar field to remain constant in some local region around it. The second scenario 
Barrow called gravitational memory because the black hole would locally preserve the value of the scalar field from 
when it formed. Barrow & Carr studied the evolution of primordial black holes for these scenarios and found 
that either case results in a significant deviation from the usual GR analysis. They considered ST theories, where 
G{(j)) oc as in the Brans-Dicke case. Since (j) increases with cosmological time, this implies that G{(j)) decreases, 
so black holes would take longer to radiate away their mass via Hawking evaporation than in the GR case. In both 
scenarios the black holes form when gravity is stronger, so the rate of evaporation is less, but in scenario (B) the 
strength never decreases and so the lifetime is even longer. 

The above scenarios are the two extremes and the reality is probably somewhere in between. Two more general 
scenarios have been proposed 0. In scenario (C) the scalar field evolves faster at the event horizon (EH) than at the 
particle horizon (PH), so 4>eh > 4’ph- Eventually the black hole must reach a stage where 4>eh = 4’ph, but this does 
not necessarily mean the scalar field is homogeneous since there could be some lag between the asymptotic and local 
increase. In scenario (D), the scalar field is evolving locally but at a slower rate than asymptotically, so (j)EH < <Pph- 
There is still some gravitational memory but not in the strict sense of scenario (B). In this scenario the gradient of 
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the scalar field is increasing but one would expect there to be some limit, due to the influx of scalar gravitational 
waves. 

Gravitational memory for black holes which are small compared to the cosmological scale (i.e. the particle horizon) 
has already been investigated by Jacobson |^. In this case, the scalar field evolution can be considered as an asymptotic 
perturbation to the Schwarzschild metric and the end-state of scenario (C) applies, with the lag being found to be 
small. This suggests that gravitational memory is virtually non-existent. However, this approximation may not apply 
for primordial black holes since these can have a size comparable to the particle horizon at formation ^ . It is still not 
clear what would happen in this case. Therefore another way of investigating gravitational memory, without assuming 
that the black hole is small, is needed. 

In GR there have been several attempts to study analytic solutions which represent black holes within a cosmological 
background. The earliest used the Einstein-Straus solution 0 , which matches a Schwarzschild interior to a Friedmann 
exterior, and this approach has also been used to study gravitational memory |^. However, in most circumstances 
it can be shown |j^ that such a matching is only possible if the scalar field is constant, which just gives the GR 
solution. Another method used the McVittie metric |^, but it has been shown |Q that this has a scalar curvature 
singularity at the event horizon. 

A more successful method uses the Tolman-Bondi metric to represent the collapse of dust to a black hole in 
an asymptotically Friedmann background. However, this only works for dust and cannot in general be applied to ST 
theories due to the derivatives of the scalar field appearing as source terms in the field equations. In this paper we 
overcome this problem by assuming that the effect of the scalar field on the spacetime is small compared to that of 
the matter. This means that we can use the usual Einstein field equations to generate the spacetime and then use 
the wave equation for the scalar field to determine its evolution. This approximation was used by Harada et al. |T^ 
to calculate the scalar gravitational radiation emitted by Oppenheimer-Snyder collapse in ST theory. Jacobson also 
uses this approximation when calculating the effect of an evolving scalar field in Schwarzschild spacetime. We note, 
however, that self-consistent numerical calculations of spherical gravitational collapse in ST theory in asymptotic 
flat spacetime, in which the effect of the scalar field on the spacetime is fully incorporated, have been considered by 
previous authors jl7[ H). 

If one makes this approximation to investigate gravitational collapse in a Tolman-Bondi spacetime, the solution is 
specified by two arbitrary functions: the energy and mass functions. To represent collapse in a Friedmann background, 
the energy function has to be negative within some radius tq and zero outside it. This results in the eventual 
gravitational collapse of all the matter within ro, while the matter outside ro expands forever as in a flat Friedmann 
model. In choosing the mass function, or equivalently the density perturbation, we adopt the “compensated” model. 
This means that the overdense region in the centre is surrounded by an underdense region outside, so that the total 
mass at infinity is unaffected. 

We solve the field equations numerically using the characteristic method. This method was first applied to an 
inhomogeneous and dynamical background spacetime by Iguchi et al. p^ . The characteristic method integrates over 
null hypersurfaces with the event horizon as a boundary. This means that one never needs to calculate anything inside 
the black hole, thereby avoiding any numerical problems associated with singularities. The output of the code shows 
the spatial and temporal variation of the scalar field. The figures produced show that the initial collapse results in 
a large gradient in the scalar field. However, as time increases, the scalar field becomes almost homogeneous. This 
suggests that, within the approximation used, gravitational memory is not possible. It remains to be seen whether 
the back reaction of the scalar field could alter this conclusion. 

In Section H we describe ST theories in more detail and derive the field equations for the approximation in which 
the effect of the scalar field on the metric is neglected. In Section HI we transform the Tolman-Bondi solution into null 
coordinates, giving the equations necessary to apply the characteristic method. In Section IV we specify the model 
giving rise to black hole formation. We present the numerical results in Section V and discuss their implications in 
Section VI. 


II. BASIC EQUATIONS 

For scalar-tensor theories of gravity with G{(j>) cx (j)~^ the field equations are 


r - -c 

^ab — I ^ ab “r ,9 
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where Tab is the usual energy-momentum t enso r, T is its trace and we have set c = 1. In Brans-Dicke theory Eq. (13) 
remains unchanged but duj/dcj) = 0 in Eq.(2.2). These equations are expressed in the Jordan frame but ST theories 
can also be expressed in a conformal frame known as the Einstein frame. The Einstein frame is related to the Jordan 
frame by the transformation 


9ab = {Go^)9ab ^ fab = (Go0)-^T,b, (2.3) 

where Gq is the present value of the gravitational “constant” as measured in solar system experiments. It is called 
the Einstein frame because it can be expressed as GR with a scalar field. However, the Jordan frame will be used 
throughout this paper. 

In GR the Tolman-Bondi solution is given by 


ds^ = -dt^ + A^{t,r)dr^ + R^{t,r)[de^ + sin^^#^], 
where r is the comoving radial coordinate and R is given by 



(2.4) 


(2.5) 


Here G{y) is a positive function given by 
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ts{r) is a constant of integration, and F{r) 


2Gom(r) with m(r) being the mass within radius r. A is given by 


and R satisfies 


A^t,r) 


l + /(0’ 




The density of the dust p is given by 


F' 
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(2.7) 
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In the above, a dot denotes dt and a prime denotes dr- 

There are two arbitrary functions in this solution: the mass function m(r) and the energy function f(r). Investigat¬ 
ing collapse to a black hole just requires the appropriate choice for these functions. In this paper the energy function 
is chosen such that 


f{r) <0 for r < To 

(2.10) 

0 

A 

0 

II 

TT 
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for some tq. This means that when a perturbation is applied to the background dust, all the matter interior to rp will 
eventually collapse to form a black hole, while the exterior matter will expand for ever as in a flat Friedman universe 
p0| . The mass function m(r) is determined by putting R = r aX t = ts{r) in Eq. (2.5). 

The key to this approximation is that the back reaction of the scalar field is neglected. It is assumed that the effect 
of the scalar field on the spacetime is small compared to that of the matter. The initial configuration used is the 
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general relativistic one with constant (j) and for simplicity Brans-Dicke theory is used. Then, to the lowest order, the 
evolution of (j) is determined by the wave equation: 


[V=VeU</'=3^rTB (2.12) 

where c]tb Ttb are determined for the Tolman-Bondi metric and the general relativistic solution. Using 
the Tolman-Bondi metric, the wave operator is given by 




{AR^y 

AR^ 


1 

I2' 


1 /n 


AR^ \ A 


(2.13) 


III. CHARACTERISTIC METHOD 

The last equation needs to be rewritten in terms of a null coordinate suitable for the characteristic method. The 
retarded time coordinate u is introduced such that u = constant is an outgoing null geodesic. In the original 
coordinates the outgoing null geodesic is given by 


dt 

dr 


= A 


(3.1) 


so we can write 


- = -A 


(3.2) 


The coordinate system is now transformed from (t,r) to (u{t,r),r{r)) using the relations 


du = — (dt — Adr ), dr = dr, 
a 


(3.3) 


where 


a = 


(3.4) 


The partial derivatives are then related by 


dt = —du, dr = -—du + df. 
a a 


(3.5) 


In this coordinate system the metric becomes 

ds^ = —a^dv? — 2aA(u,r)dudf + R^{u,f)[dO^ + sin^ Od'tp'^]. (3-6) 

To use the characteristic method it is necessary to introduce the derivative along the ingoing radial null geodesic. 
In the original coordinates the ingoing null geodesic is given by 


dt 

dr 


-A, 


which in the new coordinate system becomes 


dr a 

du 2A 


(3.7) 


(3.8) 


Therefore the derivative along the ingoing null can be obtained: 
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du^~ 

The partial derivatives dt and dr can now be rewritten in terms of d/du and df 
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The wave operator then becomes 
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ip = df{R4>). 


(3.13) 


Here ' and ' refer to the operators given in Eq.(3.1C) and Eq.( |3.11 ) respectively, 
equation for a. This is achieved by using {ii)' = (u')', which gives 


It is also necessary to obtain an 


dftt = Aa. 


(3.14) 


Applying the full Tolman-Bondi solution, the basic equations that must be solved in Brans-Dicke theory are to first 
order 
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where the upper and lower signs correspond to an expanding and collapsing phase respectively. For numerical purposes 
it is also convenient to take the parametrized form of the Tolman-Bondi solution. For / = 0, R is given by 


For / > 0, it is given by 


«.(?£) 'V, -7.M+. 


F 

R = ^(cosh77- 1), 


For / < 0, it is given by 
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Here the signs are chosen so that they correspond to the big bang universe. 
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IV. MODELS 


We choose the background primordial black hole model so that the following conditions are satisfied. (1) The big 
bang occurs at the same time everywhere, i.e. ts{r) = 0 (constant). (2) The model is asymptotically flat Friedmann 
and compensated (i.e. the overdense region in the centre is surrounded by an underdense region outside in such 
a way that the total mass at infinity is unaffected). (3) The model is free of shell-focusing or shell-crossing naked 
singularities, at least within the calculated region. (4) The central region is bound, while the asymptotic region is 
marginally bound. (5) At the initial time t = to, the condition i?' > 0 is satisfied everywhere. 

In order to satisfy the above conditions, we set t = to and choose the energy function /(r) to have the form 


/W = 



for 


for 


r <r^ 


r > Tw 


(4.1) 


where gives the cur vature radius in the central closed Friedmann region and Tw gives the scale of the overdense 
region. Equation (^^) means that the central region r < \s described by the exact closed Friedmann solution, 
which ensures that there is no shell-focusing naked singularity. More general situations in which shell-focusing is 
avoided have been discussed elsewhere . The code continually monitors for shell-crossing to check that this never 
happens within the calculated region. We then determine F{r) so that r coincides with R at the t = to spacelike 
hypersurface. 

The choice of the function f{r) requires some justification since it is nowhere exactly zero in the calculated region, 
so strictly speaking the entire region would eventually collapse to a black hole if one waited long enough. In practice, 
this does not matter because our calculated region is so large that the value of / is effectively identical to 0 for r > 

(i.e. at the outer boundary). Since the region is finite, we can therefore always make the matching to the Einstein-de 
Sitter universe outside the calculated region. It would be easy in principle to perform the calculations for a situation 
in which the matching occurs within the calculated region (i.e. with / = 0 exactly at the edge of the region). It is 
clear that this would make no qualitative difference to our conclusions but, to confirm this, we also adopt the choice 


/(?-) 



r <r^ 
r > rw 


(4.2) 


in which the compensation is satisfied explicitly within the calculated region. 

Before integrating, we have to fix the initial data for (j). One can take this to be the homogeneous cosmological 
solution given by 

We then set the initial null hypersurface as the null cone whose vertex is at {t, r) = {to, 0) and regard the cosmological 
solution as the initial data on this hypersurface. Although the value of the scalar field at the cosmological particle 
horizon must be given by this solution, the value in the perturbed region and the surrounding region may be different 
from this. To examine the sensitivity of the results to this alteration, we consider another form of the initial data 
which is different from the cosmological solution in the region t < (1 — 10) x to- We now choose 


D± = 



/1 yi' 

1 ± exp 



(4.4) 


so that we have an ingoing wave in the perturbed region, and examine the evolution of the scalar field thereafter. 
The numerical code has been checked by the following non-trivial test calculation. In the flat Friedmann universe, the 
code must reproduce the cosmological evolution Eq. (O) from the initial data. There is agreement to within 0.05% 
accuracy. 


V. RESULTS 

We denote the Hubble parameter in the Friedmann background (far from the perturbed region) as Ho at t = to- 
Recall that rc gives the amplitude of the density perturbation, while rw gives the size of the perturbed region. For 
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super-horizon scale perturbations, we cannot set the density perturbation to be very large else the overdense region 
closes up on itself and becomes disconnected from the rest of the universe j|]. Actually, the requirement R' > 0 
imposes an even stronger condition since it implies Tw ^ Tc- We have set the Brans-Dicke parameter to be oj = 5. If Tc 
is much increased, then the amplitude of the overdensity is much decreased and the resulting black hole becomes very 
small compared with the horizon scale at the formation time. If is much decreased, the overdense region becomes 
a separate closed universe. Models and parameters are summarized in Table |. The difference between models A, 
B, C and F is in the choice of the background perturbation. The difference between models A, D and E is in the 
choice of the initial cond ition for the scalar field. A change of uj only scales the variation of the scalar field from (j)Q, 
as indicated by Eq. (|4.3|). 

In Fig. 1^, the energy function /(r) and the initial density perturbation 


XU _ P{to,r) - p{to,oo) 

Pito,r) = -77-r- 

p{to,oo) 


(5.1) 


are plotted. The trajectories of outgoing null geodesics are plotted in Fig. ||. It is seen that a nearly horizon-scale 
black hole is formed for model A, while the black hole is smaller than the horizon scale for models B, C and F. This 
is because the size of the black hole is always roughly 20Hq^, whereas the time at which it forms is about 20Hq^ 
in model A and closer to 100i7(C^ in the other cases. The initial data for the scalar field are set on the initial null 
Cauchy surface (see Fig. We prepare three sets of initial data, (j)c, </>+ and </)_, and these are plotted in Fig. ||. 
We have investigated all the models listed in Table || and the results are seen in Figs. In Fig. ||, the profile of the 
scalar field is plotted for constant t. In Fig. ||, it is plotted for constant R. The reason why some of the curves come 
to an abrupt end below some value of R is that the event horizon has formed near the center and we did not calculate 
the evolution of the scalar field inside the event horizon. The reason why some of them come to an abrupt end above 
some value of R is due to the finiteness of the region in which the numerical calculation is done. 

We note that collapse ensures that the scalar field is initially concentrated in the central regions and this means 
that it rises above the asymptotic cosmological value everywhere. However, this central concentration tends to fall 
due to the underdensity surrounding the black hole. This effect, coupled with the increase of the cosmological value, 
means that the scalar field necessarily falls below the cosmological value at sufficiently large values of R, at least for 
the models under consideration. Eventually it may do so at every plotted value of R. Strictly speaking, the issue of 
gravitational memory is concerned with the process whereby the scalar field is raised to the cosmological value once 
it has fallen below it rather than with the process whereby it initially falls to the cosmological value. 

Generally, the configuration of the scalar field tends to spatially homogeneity with increasing time and the value 
of the scalar field around the black hole follows the cosmological evolution of the scalar field at each moment. To see 
this more clearly, we identify the value of the scalar field on the latest null ray within the calculation as the value on 
the black hole event horizon, 0 bHj because the null ray is very close to the event horizon for each model (at least for 
Hot < 300, see Fig. i- Both (fc and 0 bh are plotted in Fig. The end of the curve ^bh corresponds to the formation 
time of the event horizon. We can see that (()bh follows the evolution of (fc- There is a small deviation of 0 bh from 
(j)c but this can be explained by the central overdensity and the surrounding underdensity. For Hot < 100, the event 
horizon runs through the overdense region and hence the scalar field is amplified compared with (pc- Thereafter the 
event horizon runs through the underdense region and so the scalar field is reduced compared with pc- The overall 
evolution of ^bh is well described by pc- 

Although there are minor differences among the models, we can conclude that the configuration of the scalar field 
is nearly spatially homogeneous and well described by the cosmological solution pc, at least for around ten initial 
Hubble times after the formation of the event horizon. 


VI. SUMMARY 

We have calculated the evolution of the Brans-Dicke scalar field in the presence of a primordial black hole formed 
in a flat Friedmann background. We have found that the value of the scalar field at the event horizon almost always 
maintains the cosmological value. This suggests that primordial black holes “forget” the value of the gravitational 
constant at their formation epoch. In this sense, we confirm the result of Jacobsen j|], although he never carried 
out an explicit evolutionary calculation. While it is possible that some choices of the function /(r) might lead to a 
different conclusion, we would claim that our conclusion must at least hold in compensated models. For we have tried 
many sets of parameters, corresponding to models which describe a wide variety of physical situations, and this always 
seems to be the case. However, it is not clear from our analysis what would happen for non-compensated models and 
we have also avoided models with shell-crossing singularities, so we do not claim that our result is completely general. 

It should also be stressed that this result has only been demonstrated for a dust universe in which the scalar field 
does not appreciably affect the background curvature and it remains to be seen whether the same conclusion applies 
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when this assumption is dropped. As can be seen from Figs. ^ and both the radial gradient and the time derivative 
of (j) are large at early times and both of these act as source terms in the field equations. However, as long as ut is 
large, this is unlikely to stop the scalar field becoming homogeneous eventually, since it is clear that both the radial 
gradient and time derivative become small well after the initial collapse. Neglecting the back reaction should therefore 
be a reasonable assumption in this situation, so gravitational memory seems unlikely. However, this conclusion might 
not apply for a; « 1. 
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TABLE I: Parameters for models 


Models 


/ Horc HoCw initial data 
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1.25 
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1.25 

1.25 

1.25 
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FIG. 1: (a) Energy function / and (b) initial density perturbation 5 are plotted for models A-F. 

FIG. 2: Trajectories of outgoing nnll rays are plotted (a) for models A,D,E, (b) for model B, (c) for model C and (d) for model 

F. 

FIG. 3: Penrose diagram of a primordial black hole in a flat Friedmann universe. Also depicted is the null Cauchy surface on 
which the initial conditions are set for the numerical calculations. 

FIG. 4: Initial null Cauchy data sets for the scalar field: (j>c for models A,B,C and F; (j>+ for model D; and (j>- for model E. For 
clarity, the abscissa is plotted logarithmically. 

FIG. 5: Configuration of the scalar field at each moment t — const is plotted for models (a) A, (b) B, (c) C, (d) D, (e) E and 
(f) F. For comparison, the cosmological value 4>c at each moment is also plotted as a horizontal line. 

FIG. 6: Time variation of the scalar field along the world-lines of constant R is plotted for models (a) A, (b) B, (c) C, (d) D, 
(e) E and (f) F. For comparison, the cosmological evolution i^c is also plotted. 

FIG. 7: Time evolution of the scalar field on the event horizon is plotted for models (a) A, (b) B, (c) C, (d) D, (e) E and (f) 
F. For comparison, the cosmological evolution i^c is also plotted. 
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